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Abstract 

It is shown that the Lindblad equation accounts for memory ef- 
fects. That is to say, Lindblad operators can be constructed in a 
natural manner such that a memory term appears in the asymtotic 
{time — > oo) region; at the same time the expectation values depend 
on the initial state. Furthermore a procedure to extend the Lindblad 
equation to an equation of motion for an ideal Bose 'gas' of 'parti- 
cles', i.e. systems with non-trivial internal structure, is described. Initially 
in some quantum state this collection of 'particles' will asymptotically 
turn into an equilibrium ensemble whose probability distribution is 
determined by the Lindblad operators building the dissipative part of 
the equation of motion. 

A) The Lindblad generalisation ^ of Schroedinger motion to dissipative mo- 
tion hinges, apart from technical assumptions, only on very general physical 
notions: 



i) The abelian group of unitary Schroedinger motion generated by the 
Hamiltonian generalizes to a set of abelian semigroups characterized 
by a collection of operators Vj,the Lindblad operators. The semigroup 
structure - the non-existence of time-reversed motion - accounts for ab- 
sorption. 
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ii) Complete positivity, interpreted physically, means that positive mo- 
tions in Hilbert space of states (system 1) can be extended to a posi- 
tive motion in the product space resulting by adjoining a second Hilbert 
space (system2), a construction leading to entanglement of both sys- 
tems. 

It should be noted that in the derivation of the Lindblad equation concepts 
used in the derivation of master equations - for instance the decomposition 
of the space of states in a product 'system! ®' bath' - do not play any role, at 
no place. Nor is there any conceptual relation to open systems. The relation 
of master-equations for open systems and Lindblad equations has been clar- 
ified to some extent 00 and shown to be controlled by relative 
scales. 

We take up this observation and consider a system whose degrees of freedom 
interact with scale-dependent Hamiltonians and look for stationary states 
evolving from given initial states, i.e. we construct maps (p is the density 
operator of our system) 

TiV) : Q\t=o ^ Q\t=oo (1) 

and discuss their dependence on the Lindblad operators Vj which together 
with the Hamiltonian are supposed to differ for different scales - time scales, energy 
scales etc. 

We present an explicit construction of asymptotic stationary states which 
will be seen to contain memory terms. 

B)In this section we consider the case of only one Lindblad operator V and 
write the Lindblad equation of motion 

B = i[H,B] + V+BV -]^[V+V,B]+ (2) 

where H is the Hamiltonian and B an observable.^ 
Using the polar decomposition (U is a unitary operator) 

V = UVV+V (3) 



*V is assumed to be invertible; B,H are bounded operators acting in a separable Hilbert 
space for which the Lindblad equation has been proved. This fact allows us to use inter- 
changeably the notions 'operator' and 'matrix' and treat the question of dimensions - finite 
or infinite - in a rather cavalier way. 
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we rewrite (2) as ( note that the assumption of unitarity of U excludes zero- 
modes in V~^V ) 



Vv+v Vv+v 



^ [H, B] ^ + U+BU 



Vv+v ' Vv+v 



( , ^ B\/V+V + \/V+VB , ^ 
2 \VV+V W+V 



(4) 



The observation (see below) that 



is a (non-normalized) probability distribution leads us to the physically plau- 
sible assumption 

W = W{H,...) (6) 

where the dots indicate further observables commuting with H. Tracing the 
equation of motion we immediately see that the trace of the rhs of the equa- 
tion of motion vanishes identically and hence 

tr{i3W) = (7) 

or 

tr{BW) = const (8) 

{W = since W depends only on conserved quantities). Needless to say we 
tacitly assume W to be traceclass. 

In we have demonstrated the following asymptotic form for B 
i) Irreducible V 

B\t=^ = b{oo)I (9) 
where I is the unit operator in Sj. 
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ii) Reducible V, i.e. 



y^^ev; (10) 

a 

where the are matrices in orthogonal subspaces S^" of i^, yield 

5|t=^-^©6"(oo)r. (11) 



a 



In the following we consider only the irreducible case and derive 

tr{W) ^ ' 

The expectation value of the asymptotic configuration then is 

<BU.>=t.(BU<.) = ^^^||fp (13) 

for all states Qq^ i.e. the expectation value is independent of the initial state, 
no memory effects are present. We see that 

is a normalized probability distribution derived from the Lindblad operator 
Translating this result into the Schroedinger picture we derive that any 
initial state tends to Pw for t — > oo, i.e. 



riy): Qq ^ go\t=oc = Pw (15) 
for all initial states ^o- 

We now turn to the question of memory effects. To show that they can be 
incorporated we extend the Lindblad equation, without changing its formal 
content, to an equation of motion for quantum subsystems separated, e.g. 
by scales, from the system built up by these subsystems. As an example we 
could take a molecule: the subsystems are spanned by the states correspond- 
ing to the inner degrees of freedom of the atoms composing the molecule, the 
system - the molecule - is built up by the atomic states of outer shells. 
To realize this construct we endow the input matrices V and H with a direct 
product structure and, to simplify matters, choose the ansaetze 
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(16) 



H = H,^,m (17) 



V,k,H,keC (18) 

where I is the n x n unit matrix and M. is a n x n matrix sub-Hamihonian, 
identical for all sites {i, k), i.e.V and H are matrices with n x n matrix val- 
ued entries indexed by (i, A;). Our ansatz for V guarantees that, in the terms 
of our example, the Lindblad operator leaves the inner degrees of freedom 
unaffected. The observable B is written as a matrix oi nx n matrices 

B = {Biu) (19) 
The probability distribution is then 



W = [(y^v'^ ® ij =: (W,k ® l) , (20) 
V has the polar decomposition 

V=U ^V+V ={u®l^ ^V+V ® I. (21) 

Tracing the equation of motion with respect to the indices (i, k) then yields 
instead of (|^) 

fr (bW^ = mfr (WHB^ - fr (wBH^ M (22) 

which leads to 

fr {BW) = C + {wfr (wHB{t)^ - fr (wB{t)H^ dt (23) 

It has to be stated that taking the trace of this n x n matrix we obtain a 
vanishing result 

Tr (bW^ = tr^^fr (^BW^ = (24) 
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in accordance with equation (j^) .This is because we have assumed 



= [W, H] 



(25) 



and, thus 



0. 



(26) 



Following the derivation given in p| we find for the asymtotic configuration 



B\t=oo = ^(oo) ( Sik ) 
where 6(oo) is now a. n x n matrix which is read off 

b{oo) = 

1 

X 



(27) 



fr (5|i=o ^) + / (wHB{t)^ - fr (wB{t)H^ dt 

(28) 

We note the explicit appearance of a memory term. Calculating the expecta- 
tion value of -B|i=oo in some state written as 



, (0) 

^(0) = Qik 



where the g^-^ are n x n matrices, we find 



(29) 



(30) 



and observe that now the asymptotic expectation value does depend on the 
initial state in concordance with the appearance of a memory term. 



So we have seen that a simple and intuitively clear generalisation of the 
Lindblad equation to an equation for dynamical degrees of freedom of sub- 
systems leads to memory effects; the asymptotic subsystem variables given in 
equations (p3[) should be interpreted as the new dynamical subsystem vari- 
ables obtained from an asymptotic averaging procedure over those degrees 
of freedom of the total system living on 'lower' scales:the equations (^) are 
clearly seen as an elimination procedure for 'environment' variables separated 
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into a statistical average and a memory term. 



C)We now turn to the case of more than one, say N, Lindblad operators. 
Wc take N finite with the provision of eventually letting ^ oo as certain 
physical models might require. The equation of motion then reads 

B = z[H, B] + Y^ (y/BVj - \ [Vj-Vj, B] ^ (31) 

We rewrite this equation as an equation operating in a direct sum of identical 
spaces 

^N^Y.®^ (32) 



and define 



to arrive at 



Bn := BIn (33) 



Hj, := Hlj, (34) 



VN■.^^^N{VJ5JK) (35) 



Bn = i[HM, Bn] + V+BnVn - \ ^^V^, Bn\ ^ . (36) 



The polar decompositions 



Vj = Uj^Vj-Vj (37) 
lead to the polar decomposition 



^iv - Un^V+Vm (38) 

where 

Un = ( Uj5jK ) (39) 
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is unitary. Employing the same procedure as above we have 
^ B ^ 

-ON 



0^ ^Wv^ 

0^ 



(40) 



Assuming either independence of Wj on J ( does depend on J in general) 
or, alternatively, Vj positive and in analogy to (P) 

Wj:=-^ = Wj{H,...) (41) 

and taking the total trace {trN pertains to the matrix indices of the N x N 
matrices introduced above, tr^j to the operators on i^) we find 

tr^tTN (^BnWn) = (42) 



and thus 



and 



trAByWj\=^ (43) 



< tr^ {W) ^^^^ 



To illustrate this result we take Vj positive and assume the following specific 
ansatz 

Wj = {V^VY (46) 
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J = 1,2,... and put 



^ = exp ( ) (47) 



2 

The expectation value for B reaches asymptotically 

< >= }^tr„ (48) 

which is simply the expectation value of the 'particle' observable B in an ideal 
Bose 'gas' of, in the average, N 'particles' at inverse temperature (3; 'particle' 
is just a more intuitive name for the physical object dubbed 'system' up to 
now. 

This interpretation deserves further clarifications. To this end we reformulate 



equation ( ^OD as an equation in Fock space S)fi aiming at the Bose nature of 



the ideal gas to be introduced. We define 



(49) 



and 



Vf := ^©y^ 



(50) 



so that 



Wf = V+Vf = J2® {V+Vf' . (51) 
J 

In the product space Sj'^'^ we select as physically relevant states symmetric 
states which we take as superpositions of symmetric 'system' product states 
- we introduce many-'particle' boson states .The observable B is extended to 
a symmetrically operating operator 



BF = ^®{...0l0B(g)...) (52) 
J 

where B stands consecutively on all positions of the J-fold product. 
The expectation value at t = oo is, in strict analogy to (|13D and (|4^) 
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< t=oo >- [06) 

where the trace is now to be calculated with a symmetric product basis in 
j^iS)j ^Qj, g^jj J Computing this trace one encounters disconnected terms (ma- 
trix elements now pertain to S^) 

< ii\BV-^'\ii >< i2\V-^^\i2 > ■ ■ ■ < ilIV-^^Iil > (54) 



ll,...,tL 



J (55) 



All these terms sum up to the same common factor in the nominator and de- 
nominator - the connected cluster theorem - so that we reproduce (^) with 
Wj = {V^ V)'^ . We conclude that any symmetric quantum many-particle 
state composed of whatever complex quantum systems - a Bose many-particle 
state - is transported by Lindblad motion into an equilibrium ensemble with 
a probability distribution 

where 



W, = — (57) 

is an operator acting in the space of the system's states. This derivation 
is a first step towards a dissipative quantum field theory: the case of free 
fields, although we never explicitly introduced this concept. I shall return to 
the extension to more complicated cases in a forthcoming publication. 
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